Discovery of New Physics in Radiative Pion Decays? 
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Recently a strong indication for a deviation from the standard model (SM) has been obtained 
by PIBETA Collaboration. Namely, SM fails to describe the energy distribution and the branching 
ratio of the radiative decays of the positive pions at rest in the high- /low- i?e kinematic region. 

The previous experiment at ISTRA facility, testing the radiative decays of negative pions in flight 
in a wide kinematic region has alarmed about the same effect, although statistically less significant. 
The present PIBETA result indicates a deficit of the branching ratio of the radiative pion decay 
in the specified kinematic region at 8a level in comparison with SM prediction, while in the other 
kinematic regions both the branching ratios and the energy distributions are compatible with the 
V — A interaction. 

We argue that this effect can result only from a small admixture of new tensor interactions. They 
may arise due to an exchange of new spin one chiral bosons which interact anomalously with matter. 
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I. INTRODUCTION 

The standard model (SM) includes three different 
types of fundamental interactions: electromagnetic, weak 
and strong ones. The corresponding forces arise due to 
an exchange of spin-1 gauge bosons, described by four- 
vector fields Va ■ For the completeness of SM at least one 
more scalar boson is required. However, due to its weak 
coupling to the ordinary matter and its big mass, this 
particle has not been detected yet by experiments. 

Meanwhile, the existence of many other bosons has 
been predicted theoretically and tested experimentally. 
The supersymmetry, for example, suggests a rich vari- 
ety of new particles, which, however, have not been yet 
observed experimentally. 

In this letter we show that the recent PIBETA result jl| 
for the radiative pion decays (RPD) 7r+ e^i^j indicates 
that a different kind of fundamental spin-1 chiral bosons 
may be present in Nature. They have not been discussed 
intensively in literature till now. 

These particles were first mentioned as a different 
type of spin-1 bosons in ref. j2[. They appear naturally 
in the analysis of possible Yukawa couplings of spin-1 
bosons to fermion currents. In the relativistic physics, 
where two different types of spin- 1/2 fermions exist, 
namely left-handed ipL — -^(1 ~ 7^)^^ and right-handed 
■0/2 = i(l + 7^)'0, two different types of Lorentz invariant 
interactions are possible: d la gauge vector interactions 

-Cy = {gLL V'L7"V'L + gRB^ ipRi^ipn) Va (1) 

and tensor interactions 

Ct - gLR Jla'^^i^RT+p + gRL J^a^^^LT-^, (2) 

where cr"^ — |[7",7'^], T'^^ are antisymmetric tensor 
fields and g denotes dimensionless coupling constants. 

Both the vector fields Va and the tensor fields 
describe particles with spin one. However, they inter- 
act absolutely differently with matter. The gauge parti- 
cles are chirally neutral and, hence, they preserve the 



helicities of incoming and outgoing fermions. On the 
contrary, the tensor bosons carry a chiral charge that 
leads to helicity flip for the interacting fermions. More- 
over, self-interactions of the chiral bosons exist even for 
an abelian case, leading to a negative contribution into 
the /3-function p|. 

An example of the presence of the new kind of chiral 
particles in Nature is the existence of the axial-vector 
meson resonance 6i(1235), which has only anomalous 
tensor interactions ^ with quarks. The successful de- 
scription of the dynamical properties of the hadron sys- 
tems hints that the same phenomenon may have place 
at the electroweak scale as well, similarly to the gauge- 
like hadron resonances p and ai, which have served as 
prototypes for the photon and the weak bosons 5]. 

The PIBETA result confirms the anomaly already ob- 
served at the ISTRA facility |2j for the negative pion 
decays tt^ — + e^i'j in flight. It is remarkable, that these 
two experiments with absolutely different systematics re- 
vealed the same effect: the deficit in PRD yield. More- 
over, the PIBETA experiment has measured simultane- 
ously absolute total tt — > ev, — > n^e'^i', /i — > evv, 
as well as partial fi — > evi'^ branching ratios, which are 
in excellent agreement with SM predictions, with better 
than 1% accuracy The only deviation from SM is ob- 
served in the high- /low- _Ee kinem^ti^ region exactly 
where the effect has been predicted 8, 9]. The effect is 
so big that it cannot be explained within SM J^J and its 
supersymmetric extensions [T]| . 

Historically various types of local four-fermion effective 
interactions have been introduced to test the eventual 
presence of a new physics. Here we will follow also the 
phenomenological approach of the effective interactions 
because the complete theory for the interacting tensor 
particles is not constructed yet. However, we argue that 
the effective interactions should include new non-local 
momentum dependent tensor interactions, in order to de- 
scribe the PIBETA result without contradiction with the 
present experimental data. 
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II. NEW TENSOR INTERACTIONS 

First of all we define the Lorentz structure of the new 
interactions on the basis of PIBETA and ISTRA re- 
sults. Since both collaborations have observed a deficit 
of events in comparison with SM expectation, it should 
stem from destructive interference between SM and the 
new interactions in RPD. 

It is well known [T^ that the dominant contribution 
in RPD comes from the inner bremsstrahlung (IB) pro- 
cess near the edge + — rnTr/2 of the kinematically 
allowed region, which include the high-E^ /\ow-Ee kine- 
matic region, where the deficit of the events was observed. 
We assume that the new interactions interfere destruc- 
tively with this process leading to the deficit of events. 
Hence, the new interactions should have the same chiral 
structure as IB which is different from V — A interactions. 
Therefore, they do not interfere with the latter due to the 
smallness of the electron mass. This is one of the reasons 
why such type of interactions are hard to detect and why 
they have not been observed before. 

Let us consider all possible current-current interactions 
which lead to a helicity flip. There are only three pos- 
sible Lorentz structures obeying this property: scalar, 
pseudoscalar and tensor currents. The matrix element of 
the pion-photon transition for the scalar quark current 
(7|u(i|7r) is zero by kinematic reasons and does not con- 
tribute to RPD. On the other hand the contribution of 
the pseudoscalar quark current to the matrix element of 
the ordinary pion decay 
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z(0|M7'5d|7r) = ^ S.SxlO^meA (3) 

is enormously enhanced in comparison with the standard 
chirally suppressed V—A contribution, and it is severely 
constrained by the experimental data On the con- 
trary, the matrix element of the tensor quark current 
(0|u(TQ./3d|7r) is zero by kinematic reasons, and does not 
contribute directly to the pion decay, thus escaping the 
experimental constraints. 

Hence, only less constrained tensor interactions are 
possible candidates for the explanation of RPD anomaly. 
We argue that particularly non-local momentum depen- 
dent tensor interactions are responsible for the detected 
anomaly. 

In order to explain the ISTRA anomaly a new quark- 
lepton tensor interactions 

= -%/2/tG' uaapd eUapi^L + h.c, (4) 

with the effective coupling constant fr — 10~^ have been 
introduced 9], here G = GpVud- The dimensionless cou- 
pling constant /t determines the strength of the new 
tensor interactions relative to the ordinary weak interac- 
tions, governed by the Fermi coupling constant Gp. Al- 
though such tensor interactions, introduced ad hoc, can 
explain the ISTRA anomaly, it has been pointed out 0| 



that the necessary value of the coupling constant /t con- 
tradicts the ordinary pion decay tt ev. This happens 
because owing to the electromagnetic radiative correc- 
tions the pscudotcnsor quark current uaafsj^d leads to 
a generation of the pseudoscalar quark current u'j^d, to 
which pion decay is very sensitive. 

The generation of the pseudotensor term uoap^^d can- 
not be avoided for derivative free local four-fermion inter- 
actions even if we have started with a parity conserving 
quark current uaapd in eq. Owing to the identity 
uUapdR eGapVL = 0, the Lagrangian |0J effectively reads 

d-^" = -\/2/tG ucFapdL ecTapVL + h.c, (5) 

where the chiral structure shows itself in the quark cur- 
rent. 

The solution of this problem was found in ref. Q via 
introducing non-local momentum dependent tensor in- 
teractions 

C't = - V^/tG uaxpdL eaxpVL 

- V2/j.G uaxadn, — — euxpVL + h.c. 

= - V2f^G uaxad ^'^f^^ eaxpVL + h.c, (6) 

Q 

where Qq is the momentum transfer between quark and 
lepton currents, and = /t are positive dimensionless 
coupling constants. In this case the second term in the 
second row of eq. 10 is no longer equal to zero, despite 
the different chiral structures in quark and lepton cur- 
rents. And due to the identity 

UaapdL eCTapl^L = UUXadh Jy2. ^CrA/JI^L, (7) 

it can compensate the opposite chiral quark structure of 
the first term. Then the terms with the pscudotcnsor 
quark currents uaapj^d cancel out in eq. ((^Jl, and the 
tensor current uaapd does not contribute to pseudoscalar 
pion decay because of parity conservation in electromag- 
netic interactions. 

The two different terms in the effective Lagrangian © 
come from exchanges of new spin-1 bosons T^^ and U^^ 
with opposite chiral charges, which are necessary to avoid 
a chiral anomaly j^. The peculiar Lorentz structure of 
the new tensor interactions reflects the Lorentz structure 
of the propagators for the chiral bosons. This structure 
can be obtained, for example, from the one-loop radiative 
corrections to the self-energy of the tensor fields using the 
interactions (0) . This follows from the fact that in case of 
dimensionless coupling constants the theory is formally 
renormalizable and the radiative corrections should re- 
produce the Lorentz structure of the kinetic terms in 
the bare initial Lagrangian for the tensor fields. 

In general, the coupling constants /t and /-^ can be 
different, however, we assume their equality to avoid the 
experimental constraint from the ordinary pion decay. 
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In the following we will keep different notations for the 
coupling constants in order to compare the effects from 
the two different Lagrangians Q and ©. 



III. RADIATIVE PION DECAY 

The most general matrix element of RPD tt^ —^ e^v^ 
reads 



M = MiB + MsD + Mt + Mt: 
where besides SM matrix elements for IB process 

_ pq 



MiB = i—^J^rrieea e 



kq 



(8) 



'^L (9) 



and for the structure-dependent (SD) radiation 
.V2eG 



MsD 



£a { FA[{pq)gat3 -Paqfs] 



the new tensor contributions 

Mt = -V2eGFT Saqp ecTapvL 

and 



M't 



-V2eGFT {qc 



^ Q\Ql3 - 

£\ - q\£a) —7yr~ ^^apVL, 



(10) 

(11) 

(12) 



are present. Here is the photon polarization vector; 
p, fc, and q are pion, electron and photon momenta, cor- 
respondingly. 

The first term in eq. © Mi b describes a gauge invari- 
ant QED process of the photon radiation from the ex- 
ternal charged particles, which is free from the effects of 
the strong interactions. It contains only one well-known 
phcnomcnological parameter: the pion decay constant 
= (130.7 ±0.4) MeV. 

The second term Msd corresponds to the photon 
emission from hadronic intermediate states, governed 
completely by the strong interactions physics. It is 
parametrized by the two form factors Fy and Fa of the 
7r-7 matrix elements for the vector quark current 



{l{q)\uia.d\lT (p)) = 80 FyeappaPp qa 

and for the axial-vector quark current 



(13) 



{'y{q)\ujai d\Tr (p)) = i — Ef^FA [{pq)gaf3 - qapp] 



(14) 



Assuming CVC hypothesis, the vector form factor Fy 
is directly related to the tt" — > 77 amplitude 'iS] and can 
be extracted from the experimental width of this decay 



Fv 



l2T{-K° 77) 



0.0262 ±0.0009. 



(15) 



This value is in a fair agreement with the calculations 
in the relativistic quark model p/7| and with the lead- 
ing order calculations of the chiral perturbation theory 
(CHPT) [13 



Fv 



1 

47r2 F^ 



0.0270. 



(16) 



The value of the axial form factor Fa is model depen- 
dent and its determination is a matter of experimental 
measurements. The ratio of the axial to the vector form 
factors 7 = FaI Fy has been measured in the previous 
experiments [l9l | in kinematic regions where the contribu- 
tion of the new tensor terms is not essential. The average 
value 7 = 0.448±0.062 at fixed Fy = 0.0259±0.0005j20| 
is also in agreement with the calculations in CHPT [2l|. 

The matrix elements Mt and M^ follow from the new 
interactions between quark and lepton tensor currents 
©. The matrix element for the quark tensor current 

{l{q)\u<7apl^d\-K-{p)) = -^Pt iqa£p ^ qfiSa) ■ (17) 

can be calculated [m applying the QCD sum rules tech- 
niques and the PCAC hypothesis. So 



(18) 



3 /tt 



is expressed through the magnetic susceptibility [22, |2j{ 
X = — (5.7± 0.6) GeV~^ of the quark condensate and 
its vacuum expectation value (OlggjO) w —(0.24 GeV)'^. 
Then the tensor form factors in eqs. ^\\\ and (|12|l read 
Ft = -{fT + f^)F° and F^ = -2/^^°. 

In general, all form factors depend on the square of 
momentum transfer to the lepton pair — (j> — qY ■ 
However, these dependences are weak and, hence, the 
form factors can be assumed as constants. 

The differential decay width of RPD 



d^r 

dxdy 



27r [l-ry 



Pi.x,y) 



(19) 



can be expressed through the ordinary pion decay width 
^■K^ev, where the kinematic variables x — 2pq/m'^, y = 
2pk/rn^ and the ratio r — (me/mjr)^ ~ 1.34 x 10~^ are 
introduced. The Dalitz plot distribution is defined by the 
density 

p{x,y) = piB{x,y) + psD{x,y) + piBSD{x,y) 

+ PT{x,y) + psDT{x,y) + piBT{x,y), (20) 



where 

PIB 
PSD 
PIBSD 

Pt 

PSDT 
PIBT 



IB{x,y), 

[{\+^fSD+{x,y) + {\ ~ SD' (x,y)\ , 
2a^[{l + ^)G+{x,y) + {l--f)G-{x,y)\ , 
T{x,y), 

2a^V^ [(1 + 7) J+(x, y) + {l~ i)J- (x, y)] , 
2al{x,y). (21) 
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The explicit forms of the functions IB{x,y), SD'^{x,y), 
G^{x,y), T{x,y), J^{x,y) and I{x,y) are given in the 
Appendix. The constant 



■Fv = 



3.945 



(22) 



defines the strength of IB contribution relative to other 
contributions. 



IV. DISCUSSION AND CONCLUSIONS 

Based on the previous consideration, we discuss now 
the experimental data and their interpretation. To inves- 
tigate the most interesting part of RPD, namely the SD 
radiation, and to extract 7, all previous experiments 
have been fulfilled in a restricted kinematic region com- 
patible with region A of PIBETA experiment (Fig. ^), 
which is an intersection of regions B and C. 
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FIG. 1: (a): SM distribution (ISTRA region is marked out by 
dashed lines); (b,c,d): isocurves (with 10% step) correspong- 
ing to relative strength of tensor contribution for different 
values of Ft and F^ (the dashed curve corresponds to zero). 

In the Fig. ^ the isocurves for pure SM contributions 
(Ft=F!j. = 0) are shown. IB contribution dominates near 
the edge x -\- y — \ oi the kinematically allowed region 
and SD+ reaches its maximum near the point (2/3, 1), 
while their interference and SD~ contribution are small. 

Almost the whole kinematically allowed region 
(Fig. ^i) has been investigated first at ISTRA facility . 
A large deficit of events (33 ± 10)% in comparison with 
SM prediction has been observed. Even with a poor 



statistics they were able to establish the kinematic region, 
where a lack of events occurs. This region corresponds 
to the bottom part of region B of PIBETA experiment, 
where only the IB and the small SD^ contributions are 
expected. 

The introduction of the tensor matrix element Mt 
with Ft — —0.01 "2^ explains the lack of events and 
leads to a considerable improvement of the y^. However, 
such a big tensor form factor contradicts [Tjl the present 
experimental data. Moreover, in order to explain the 
(19.1 ±2.5)% deficit of events in region B of PIBETA ex- 
periment, using the same form (|ll|l for the tensor matrix 
element, an order of magnitude smaller value of the ten- 
sor form factor Ft ~ —0.0016(3) is required [2^. These 
different Ft values indicate an inadequate description of 
the new interaction. 

To analyse the problem further, we compare the rel- 
ative strength of the tensor contribution with respect 
to SM contribution for two different values of the ten- 
sor form factor Ft (Fig. ^ and^). The main differ- 
ence between these two plots is the following: At the 
biggest value of \Ft\ both negative and positive contri- 
butions are present depending on the region, while for 
\Pt\ ^ « 0.0069 the tensor matrix element 

leads to the only negative contribution into the whole 
kinematically allowed region. 

As one can see from the Fig.^ the tensor contribution 
with the biggest value \Ft\ affects significantly regions A 
and C, besides region B. The latter is in contradiction 
with PIBETA results. The lowest absolute value of 
allows to evade this problem but leads to another latent 
difficulty of 7 determination. 

Indeed PIBETA claims two different values for the ra- 
tios of axial to vector form factors 7 = 0.443 ± 0.015 and 
7 ^ 0.480 ± 0.016. The first one corresponds to the SM 
fits l] to the entire data set, while the second - to region 
A data only. However, the Sti deviation in the branching 
ratio from SM prediction in region B requires an intro- 
duction of the new tensor contributions for proper 7 de- 
termination. Then taking into account the correction for 
the destructive interference, the experimental branching 
ratio in region A should be increased. This increases the 
7 value, as follows from the top panel of Fig. 4 in ref. 
However, if we believe in CHPT calculations 7 should 
be decreased in order to approach the lowest value. In 
other words the fits made with Flp = are inappropriate. 

The case Ft = F^ ~ —0.01 allows to describe both 
ISTRA and PIBETA anomalies as well, without contra- 
diction to the experimental data. The corresponding ten- 
sor contribution has a slight slope near its zero value 
(Fig. ^) in regions A and C which is in accordance 
with PIBETA results. Moreover, contrary to PIBETA 
fit, mainly the positive contribution in region A leads to 
the right direction for 7 correction, corresponding to de- 
creasing its value. While the tensor contribution can lead 
up to 40% deficit in region B. 

The nonzero form factors Ft and F^ indicate an exis- 
tence of the quark- lepton tensor interactions © with the 



coupling constant fx = — 0.013. Although such type 
of interactions can be generated through radiative correc- 
tions, it is impossible to get the tensor coupling constant 
to be larger than 10"^- 10"*^ in SM and 10"''- lO'^ in its 
SUSY extensions JJJ. Moreover, the particular form © 
cannot arise as a result of Fierz transformations from a 
leptoquark exchange |26| as well. The only natural source 
to produce the effective interaction 10 is the exchange of 
the new chiral bosons interacting anomalously with 
matter. 
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yi 



[{x - yi){l -x)- rx] , 
- 2;) - (1 - r)x] , 
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It 



1-x 



X - 2yi 



1-x 



x~yi 



l-x 



ItVi + 2(7T - lT)yi 



x X — yi 



,(23) 



Appendix 



The analytical expressions for the functions IB(x,y), 
SD^{x,y), G^{x,y), T{x,y), J^{x,y) and I{x,y) read 



IB{x,y) 



a;2(x - yi) 



1 - r 



+ 2(1 - x) - 



2rx 



x-yi 



SD+{x,y) = {x-yi)[{x-yi){\-x)-rx\, 
SD~{x,y) = yi [yi{l - x) + rx] , 



where yi = 1 — y + r, jt ~ Ft / Fy and 7^ = F!p/F\/ 
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